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Abstract 

The extraction of the nucleon's strangeness axial charge, As, from inclu- 
sive, quasielastic neutral current neutrino cross sections is studied within the 
framework of the plane-wave impulse approximation. We find that the value 
of As can depend significantly on the choice of nuclear model used in analyz- 
ing the quasielastic cross section. This model-dependence may be reduced by 
one order of magnitude when As is extracted from the ratio of total proton 
to neutron yields. We apply this analysis to the interpretation of low-energy 
neutrino cross sections and arrive at a nuclear theory uncertainty of ±0.03 
on the value of As expected to be determined from the ratio of proton and 
neutron yields measured by the LSND collaboration. This error compares 
favorably with estimates of the SU(3)-breaking uncertainty in the value of As 
extracted from inclusive, polarized deep-inelastic structure function measure- 
ments. We also point out several general features of the quasielastic neutral 
current neutrino cross section and compare them with the analogous features 
in inclusive, quasielastic electron scattering. 
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I. INTRODUCTION 



Inclusive, quasielastic (QE) weak neutral-current (NC) reactions have received consider- 
able attention recently as a means of probing the strange quark content of the nucleon [|l|. 
In particular, analyses of the Brookhaven experiment E734 [0-H] have generated bounds on 



the nucleon's strangeness axial-vector matrix element which are essentially consistent |T0 
with the values for As, the strange-quark contribution to the nucleon spin determined from 
polarized deep inelastic scattering 0-0. Less stringent bounds on the strangeness vector 
current form factors have also been extracted from the Brookhaven data [0,3 . In a similar 
vein, one expects a determination of the ratio of proton and neutron yields in the LSND 
experiment at Los Alamos ]ll| to produce even more stringent limits on some of these form 
factors |I2| , pr5[ . The results from these studies should complement the results from several 



parity-violating (PV) elastic and QE electron scattering measurements currently underway 
at MIT-Bates |1|,|T§ and planned for both CEBAF [0-|T§ and Mainz Indeed, this 



program of semileptonic NC scattering measurements affords one with a unique, low-energy 
window on the non- valence quark structure of the nucleon [jl|. 

Since the neutrino NC reactions of interest require the detection of a final-state nucleon 
knocked-out of a nuclear target, a proper interpretation of the results in terms of single- 
nucleon physics requires that one have a sufficiently reliable understanding of the nuclear, 
many-body impact on the neutrino cross sections. Thus far, nuclear calculations have been 
performed using a relativistic Fermi Gas (RFC) model and a continuum RPA approach 
||r3| . In the present paper, we place these studies in the context of a more general framework 
by pointing out several features of QE neutral current scattering not previously realized 
in the literature. Specifically, we note the complementarity of inclusive, QE NC electron 
scattering, in which the outgoing electron is detected, and QE NC neutrino scattering, 
in which the detected particle is a nucleon. These two processes - to which we refer as 
t-inclusive and w-inclusive scattering, respectively - explore different regions of the two- 
parameter missing energy (S) and missing momentum (p) space. Consequently, in the plane- 
wave impulse approximation (PWIA), these two types of QE NC scattering may display 
different sensitivities to the many-body physics which enters the one-body spectral function, 
Sip,£). 

We use this framework to arrive at the first (to our knowledge) published estimate of the 
nuclear theory uncertainty associated with the extraction of As from the BNL and LSND 
measurements. Our approach in doing so is the following. By varying the nuclear model 
used in analyzing the QE cross section, we change S{p, £) and, consequently, obtain different 
extracted values for As. To be as conservative as possible, we choose two simple, tractable 
models lying near the extremes of the spectrum of reasonable nuclear models. Specifically, 
we employ the RFC model and a hybrid model (HM) involving harmonic oscillator shell 
model wave functions for the bound nucleons and plane waves for the continuum states. 
The former represents the "maximally unconfined" extreme, since it employs plane waves 
for bound and continuum single-particle states; additionally, as discussed below, only the 
on-shell electroweak current matrix elements occur [^]. The latter is "overconfined" in the 
sense that a harmonic oscillator basis is used, implying that the long-range behavior of the 
bound single-particle wave functions is gaussian rather than exponential as should be the 
case with finite potentials. Our approach in the present work is to model only these extremes 
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to explore the "worst-case scenario" for extracting As from QE neutrino scattering. 

In fact, these two models reproduce rather well the experimental QE response for inclu- 
sive electron scattering, especially for integrated quantities such as the Coulomb sum rule 



2T| or the responses discussed in the present work, despite the significant differences found in 



the behavior of the respective spectral functions in {£,p) space. We expect that the spread 
in extracted values of As obtained using more sophisticated treatments of the response, 
such as those which include the effects of correlations and more realistic single-particle wave 
functions, will be reasonably characterized by the difference between the RFG and HM val- 
ues. We take as the nuclear theory error, (5nuc(^'5), the difference between As(RFG) and 
As(HM). We find that at the kinematics of the LSND experiment, ^nucl^-^) ~ ±0.25 when 
the individual proton or neutron knockout cross sections are used. This error is roughly 
as large in magnitude as the average value for As determined from the deep-inelastic mea- 
surements. If, however, one considers the ratio of proton to neutron yields rather than 
the individual cross sections, as has been proposed for the interpretation of the LSND data 



T^,|T3|, we find that the magnitude of ^nucl^-s) is reduced by more than an order of magni- 
tude. In this case, the nuclear theory uncertainty is significantly smaller than estimates of 
the theoretical SU(3)-breaking uncertainty in the deep inelastic values for As p2| , |23[| . From 



this standpoint, our results complement those of Ref. |T^ which analyzed the impact of 
final-state interactions (FSI) on the extraction of As from QE neutrino cross sections and 
which found that use of i?^, as compared with the individual cross sections, significantly 
reduces one's sensitivity to distortion in the outgoing nucleon's wave function. 

In the remainder of this paper, we discuss these features in more detail. A reader 
primarily interested in the application of our analysis to the extraction of As is encouraged 
to read Sections and 0, along with Eqs. ( ^^ - |35| ) of Sec. |T|. The reader interested as 
well in the general features of u- and t-channel QE cross sections is directed to Sec. 0, 
where we discuss the difference between QE neutral current neutrino and electron scattering 
in general terms; Sec. |T|, in which we consider the implications of this difference for the 
(PWIA) analysis of QE scattering; and Sec. where we specify further to the RFG. Section 
|VII| summarizes our conclusions and is followed by an Appendix. 



II. INCLUSIVE t- AND it-CHANNEL SCATTERING 

The leading order exclusive QE cross section is generated by the Feynman amplitude 
associated with the diagram of Fig. |1|. Here, a lepton i scatters off an A-hodj nucleus to 
a final lepton state i' via the exchange of a vector boson V. In the scattering, a nucleon 
N is knocked out leaving behind an [A — l)-body daughter nucleus generally in an excited 
state. We let K'^ = (e, k) and K' = (e', k') denote the initial and final lepton momenta, 
respectively, = {En,pn) the four-momentum of the ejected nucleon and —p the three- 
momentum of the recoiling daughter nucleus. 

Following EM we define the missing energy S as 



S^^p^ + Mf_, - ^p^ + Ml_, , (1) 

where M^-i and M\_^ are the masses of the recoiling nucleus in its ground and excited 
states respectively. Thus the missing energy used in this work actually corresponds to the 
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excitation energy of the residual nucleus in a frame where it is moving with a momentum 
—p. The conditions for four-momentum conservation in the laboratory frame (target nucleus 
at rest) read 



e + MA = e +En + £ + ^P^ + Ml_^ (2a) 
k = k' + pn — p . (2b) 

If we further introduce the nucleon kinetic energy Tn = — mjv, the nuclear recoil energy 



Ta-i = yp"^ + ^1-1 — Ma-1 and the positive nuclear separation energy Es = Ma-i + ttin 
Ma, then one may rewrite Eq. (2B) as 



e = e' + Tn + S + Ta-i + Es . (3) 

At this point, we refer to Eqs. (^) and (|^) along with Fig. |I|to describe the difference be- 
tween QE (e, e')N and QE {u, N)v' kinematics. In the former case, the initial and final lepton 
energies and three-momenta are fixed experimentally. Thus, one may specify a given value of 
the energy transfer uj = e — e' and three-momentum transfer q = k — k' hy properly selecting 
the beam energy and momentum and the lepton detector settings. The experimentally-fixed 
Lorentz invariant for this process is = uj^ — < Q and we correspondingly refer to this 
process as "t-inclusive" scattering. Further, by combining Eqs. (|2b|) and (^, we may obtain 
the following relation between the missing energy and recoil momentum: 



£{p) =uo- Ta-i - Es 



(4) 



where the quantity in the square brackets is just T/v. Note that Ta-i also carries a depen- 
dence on p"^. However, for typical targets in QE scattering (e. §.,^"^0), Ta-i is much smaller 
than the other energies involved in the problem so that one may usually neglect it without 
introducing any significant error. Eq. (^ actually defines a continuous family of curves, 
parameterized by the angle between p and q. This family is accordingly bounded by the 
curve we denote by corresponding to cos{p ■ q) = —1 {p and q anti-parallel) and the 
curve S"^, for which cos(p ■ q) = 1 {p and q parallel). 

To get the t-inclusive cross section, one must integrate over the three-momenta of the 
undetected particles (daughter nucleus and outgoing nucleon). As we show below, this 
corresponds to integrating over the region in the {S,p) space lying between the two curves 

and . In fact two situations may actually occur, depending upon the sign of the 
vertical intercept of these curves: 



= uj - Ta-i - Es - {\lq^ + m% - rriN) ■ (5) 

For If < 0, the integration region, which we denote T>, will be bounded by S~{p) and the 
p-axis (Fig. ^a), while, for It > 0, V is bounded by £~{p), £^{p) and the p-axis (Fig. ||b). 
Note that in the figures, as an orientation, we have taken typical values for q (500 MeV/c), 
UJ (100 and 170 MeV, respectively) and Es = 8 MeV. Furthermore, we have included only 
the portions of the curves existing in the upper right quadrant, since both £ and p are by 
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definition non-negative. For this region, the maximum value of the missing energy occurs 

for 



Ma-i 

P = Pmax = \ (6) 

Ma-1 + rriN 



and has the value 

8 iPmax^ = 

= uo — Es - 



^ + {Ma-i + rriNf - {Ma-i + m^) 



^uj-Es . (7) 

For typical QE kinematics, one has Emax ~ 100 MeV. 

Parenthetically, we note that the region V depends on the independent kinematic vari- 
ables for the inclusive process. For example, when the cross section for an outgoing electron 
of a specified energy and scattering angle is considered, one has T) = T>[e',6e)- When the 
energy- or angle-integrated cross section is of interest, the region V consists of the union of 
all regions T>{e', Of.) allowed by the kinematics. 

In the case of QE NC neutrino scattering, the final lepton is undetected and the four- 
momentum of the outgoing nucleon is specified Q. Thus, the experimentally-fixed variables 
become 

Mo = e — E]si (8a) 
u = k — pn . (8b) 

The corresponding Lorentz invariant in this case is f/^ = — and we refer to this 
process as " u- inclusive" scattering. From the conservation relations (Eqs. (^^ and (pbD), 
one has 

uq = £ + t + Ta-1 + Es - ttin (9a) 
u = k' -p . (9b) 

Since the neutrino is massless, one has that e' = \k'\ = \u + p\, so that Eq. (^B|) may be 
re-written as 

S{p)=Uo-Es + mN-TA-i-\u + p\ . (10) 

Eq. ( p!OD is the u-channel analog of the relation in Eq. (|[). Neglecting the small p^- dependent 
recoil kinetic energy and defining a quantity 6 as 

6 = Uo - Es + rriN , (11) 

we obtain the following boundaries for the region V over which one should integrate to get 
the M-inclusive cross section 



In practice, detectors include the full solid angle for the outgoing nucleon. 
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£^ = 5 — {u + p), u and p parallel (12a) 

g- — ^ I . _ { £^ = 5 — {u — p) u and p anti-parallel, u > p (I2h) 
^ \ = 6 — {p — u) u and p anti-parallel, u < p 

We note that the reason for the appearance of two cases for £^ is that the undetected 
particle in this instance is massless. Consequently, e' - and therefore £ - depends linearly 
on + In the t-channel case, the un-detected outgoing particle is massive, rendering the 
dependence of £^ on |qr+p| quadratic and the relative magnitudes of q and p inconsequential. 

To assess the extent and the shape of the region T) it helps to observe that in the it- 
inclusive scattering the maximum of the curve £^ occurs for 



where £ {p) assumes the value 



P = Pmax = U (13) 



£ (Pmax) = S . (14) 



The relations in Eqs. ( [T3| ) and (|T^) are the w-channel versions of Eqs. (|^) and (|^). In both 
cases, the location of the maximum depends linearly on the magnitude of the independent 
vector quantity {q or u) while the height of the maximum is linear in the independent 
time-like quantity {u or uo). 

For future reference, we specify the two values of the momentum {p~ and p~) where £~ 
vanishes. They are given by 

p^ = u — 6 and p^ = u + 6 , (15) 

p~ being positive definite whereas p^ can be either positive or negative. On the other hand 
£^{p) vanishes for 

p+ = S-u , (16) 

which can be positive or negative. 

Again in complete analogy with the t-channel the curves £~{p) and £^{p) intercept each 
other at p = 0, where their common value is 

Iu = £'{0) = £^{0) = 6-u . (17) 

Let us denote by V{On) the allowed region in the {£,p) plane at fixed 6'Ar, the angle between 
k and Pat.Q Two situations can then occur corresponding to two different shapes for this 
region: Either lu is positive or is negative. In the latter case p^ is positive and V{6n) is 
a triangle, lying of course in the physical quadrant of the {£,p) plane. The area of the 
triangle is S"^; hence, it is fixed once the moduli of the momenta of the incoming neutrino 



^As in the t-channel case, the integration region depends on the two independent kinematic 
variables. In the present context, it is useful to choose one of them to be 9n and to suppress the 
dependence of T> on the other. 
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and outgoing nucleon are given. On the other hand the position of the triangle in the {S,p) 
plane depends on 9]y, as well as on A; = |A;| and = \p]\f\ (indeed, according to Eq. (|T3|), 
the upper vertex of the triangle is fixed by u, whose value also depends on these variables 
(see Eq. (^)). 

This situation is illustrated in Fig. |^a for kinematics typical of the LSND experiment, 
namely by choosing e = 200 MeV and Tn = 60 MeV {pn = 341 MeV/c). We thus see that 
when Pn and k are antiparallel [6^ = 180°), /„ is indeed negative and V is given by the 
triangle on the right-hand side of the figure. 

As 6n is varied from antiparallel (180°) to parallel (0°) while keeping k and p^ fixed the 
value of u decreases and the triangle moves continuously in the leftward direction. When 
is positive, then T>{9n) becomes quadrangular, as displayed in Fig. |b for the case in which 
k and Pn are parallel and for kinematics typical of the BNL experiment. The area of the 
quadrangle is given by 2u5 — v?. 

For an experimental situation in which nucleons are detected over the full 47r solid angle 
(as in LSND), the global integration region V will be given by taking the union of sub- 
regions V{6i\f) for all 6^. For example, as illustrated in the two panels of Fig. ^ for cases 
with /„ < and > 0, respectively, the regions V are defined by the lines whose vertices 
are labeled ABCD {1^ < 0: quadrangle) and ABODE {1^ > 0: pentagon), respectively. 



III. PLANE- WAVE IMPULSE APPROXIMATION 

Much of the interest in the inclusive (e, e')N and {u, N)v' scattering reactions stems from 
an interpretation of the cross sections employing the PWIA. Invoking the PWIA corresponds 
to modeling the shaded vertex in Fig. |l] as illustrated in Fig. ^. Here, one assumes that only 
one nucleon participates in the scattering by absorbing the virtual vector boson V , the 
remainder of the nucleus acting just as a spectator. Three- momentum conservation in the 
laboratory system requires that the struck nucleon has momentum p since the initial nucleus 
is at rest and the daughter nucleus recoils with momentum Pa-i = —P- Before absorbing 
the V , the struck nucleon has energy E and in general does not lie on the mass-shell. In 
the limit that final-state interactions are neglected, as is indeed the case in the PWIA, the 
outgoing nucleon is assumed to be on the mass-shell with an energy E^q = \Jp% + mj^. 
Under these assumptions, the differential cross section can be written in terms of kinematic 
and phase space factors, the square of the invariant amplitude for scattering of the incident 
lepton from a single nucleon — half-off-shell, since the struck nucleon is in general off-shell 
p5| — and a function S{p, S), referred to as the spectral function, which carries information 
on the probability of finding a nucleon inside the nucleus with momentum p and energy 

E = Ma-^P^ + Ml_, - S , (18) 

It reads: 

_ 1 1-^2 d^k' d^pN d^PA-i 
^ ~ AkMA (27r)32e' (27r)32^^ {2iif2E\_^ 

X (27r)^5W(ir + Pj,-K'-Pm- Pa-i) (19) 
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where P4 = (M4, 0) is the four-momentum of the target nucleus, Pa-i = (-E'a-i, — p) is 
the four-momentum of the daughter nucleus, A4 is the invariant lepton-nucleus scattering 
amplitude, and where the bar over A4 denotes the appropriate average over initial spins and 
sum over final spins. Performing the integral over Pa~i gives 

1 -nr-rrs d^k' d^p^ 27r 



AkMA (27r)32e' (27r)32Ejv 2E^_i 
X 6{e + MA-e - En- E^^^) . (20) 



The general form for the square of the invariant amplitude is 



\M\=g'Dym'L,,W^'^ , (21) 

where g is the strength of the fermion-vector boson coupling, Dy{Q'^) = (Q^ — My + ie)~^, 
My is the electroweak vector boson mass, L^,^ is the usual leptonic tensor appearing in 
semi-leptonic scattering and Wfj^i, is the corresponding hadronic tensor. They are given by 
the following well-known expressions: 

Lf,u = -^Tr [u{k', s')T^u{k, s)u{k', s')T^u{k, s)] (22a) 
o 

W,. = {A\Jl\f){f\l\A) , (22b) 

where 1^4) is the initial target nucleus, |/) is the final hadronic state, is a Lorentz structure 
associated with the lepton-vector boson vertex, Jx is the hadronic current operator, and 
where the appropriate average over initial and sum over final hadronic states is understood 
in Eqs. (pOD-(p2b|). Furthermore, in Eq. (|22b|) the spinors are normalized according to the 



Bjorken and Drell conventions [26| and the weak vertices are given by the usual expressions. 



In the PWIA framework only the one-body component of the nucleonic current, namely 

xuN{k', s')r^u{k, s)d\k', s')d{k, s) , (23) 

is kept. In the above a"^(/c,s) and d{k,s) are, respectively, the operators that create and 
annihilate bound nucleons {i.e., in momentum space, having Fourier component k and spin 
projection s) and F^ is the Lorentz structure associated with the single-nucleon current 
matrix element |2^. Normalizing the states according to Ref. p6| . 



{p\q) = i27rf2Ep 5ip - u) , (24) 
and substituting the expression in Eq. ( p3D into Eq. ( p2b| ) yields 

w^.- = E (^) - - mp, s)\A) 

s,s' \ V / 

X w^^{p,Pn) , (25) 

where 
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u{p, s)Tf,u{pN, s')u{pn, s')T^u{p, s) (26) 

is the so-called single-nucleon tensor and the ket |yl — 1) represents the daughter nucleus 
in either its ground state or one of its excited states. The replacement of Ea-i with 7i = 
H — -E^li inside the ^-function appearing in Eq. (^) allows one to perform a sum over all 
daughter-nucleus states by making use of the closure relation. One then obtains 

^ s,s' 

x{A\d\p,s)6{e + MA-e' -EN~E'^_,-n)d{p,s)\A)- "^'^^ 



(27r)32e' {2tt)^2En 



(27) 



Noticing that E = En — e + e' and defining the chemical potential /j, = Ma — E^_^ (not to 
be confused with a Lorentz index), one has 



where the spectral function is 

"^^^'^^^ (2^(2^; \2Ma, 

X {A\d\p,s)6{E + n-fi)d{p,s)\A) . (29) 

A further elementary elaboration of Eq. (^) leads to the exclusive cross section 

rfV 111 



dQ'dk'dQNdEN (27r)2 2k 2E 



xL,,w^%p,PN)Sip,£)^^ . (30) 

From Eq. ( ^8|) or (|30|) one obtains the t-channel (w-channel) inclusive cross section by in- 
tegrating over the undetected nucleon (lepton) momentum. In so-doing, it is convenient to 
convert to integrals over the variables p and S. To this end, we first make use of the fact 
that Ptv = Q + P and k' = u + p, so that for fixed q one has d^p^ = d^p and for fixed u 
one has d^k' = d^p. Second, we write d^p = d(f)d cos Op"^ dp, where {9,(f)) are defined with 
respect to q in the case of t-channel scattering and with respect to u in the case of w-channel 
scattering. Third, we make use of the {S,p) relations in Eqs. (|^) and (^ to transform from 
dcos9 to d£, getting 

dcos9= I (^^/P^)^^ ' ^-channel 
1 {e'/pu)d£ , M-channel 

respectively. Finally, we obtain for the cross section 
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da = (27r) J pdp J dSS{p,S) g'DviQ'y 



X L^^w{p,p^r^ , (32) 



P 



where 



dQf _ 1 d^k' 



p q (27r)32e' 
dQf 1 d^pN 



u (27r)32E 



N 



t-channel (33a) 
u-channel , (33b) 



"D denoting the allowed region of integration in the {S,p) plane (for the case of electrons, 
we have assumed the extreme relativistic limit {k = e)). 
From the above, the explicit expressions 

"^'^ ^^.--yDviQr (34) 



and 



dQ'de' 277 32 e g 

X J^pdp j —S{p,£)L^^w^'''{p,pn 



d'^a 1 1 pn 1 . 
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d^lj\fdEj\f 27r 32 e u 
X J^pdp J ^S{p, S)L^,w'^^{p, pN)Dv{Q' f (35) 

for the t- and w-inclusive cross sections, respectively, follow. Note that the azimuthal in- 
tegration has already been carried out, introducing an azimuthally averaged single- nucleon 
tensor (or, equivalently, single-nucleon cross section) as done, for example, in Refs. [p4| , pO| . 

The formula given in Eq. (0) constitutes the conventional starting point for the analysis 
of inclusive QE (e, e') scattering in the P WI A .0 Previously reported treatments of QE (z/, N) 
scattering, however, have not made use of this framework |T^J^. The advantage of writing 
dalu, N) in the form of Eq. ( ^2]) is two-fold: First, it makes explicit the dependence of the 
QE cross section on the experimentally-fixed kinematic variables (mq, u) via the specification 
of the integration region V and the appearance of u in the integration measure 1/p. Second, 
it makes the role of the one-body nuclear spectral function explicit, thereby making the 
nuclear model-dependence of the cross section more transparent. 

When the single-nucleon tensor w^^{p,pi^) refers to an on-shell nucleon, the expressions 
in Eqs. ( ^21) and (|3^) carry a dependence on the free-nucleon form factors. For general 
neutral current scattering, three such form factors contribute: Gi(Q^), where i = E,M, or 
A denotes the Sachs electric and magnetic form factors and the axial-vector form factor, 
respectively M. In the case of t-channel scattering, one is then able to employ Eq. (|33) to 



^Typically, one sees integrations over E rather than £ in the literature. 
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extract information on the single- nucleon form factors at a single value of . Since t is 
fixed experimentally in this case, one may factor the form factors out of the integral over 
0), leaving only an integral over the spectral function and various kinematic factors 
resulting from the contraction of L^^, and w^" . To the extent that one's choice of S{p,S) 
is realistic and that the PWIA is valid, one obtains a more or less reliable determination 
of the Gi{Q'^). An important feature of t-channel scattering is that g, uj and the electron 
scattering angle 6e can all be fixed in an inclusive measurement, allowing (in the plane- wave 
Born approximation) the various longitudinal (L), transverse (T), . . . hadronic responses to 
be separated before attempting to determine the Gi{Q'^), to be contrasted with the situation 
discussed below for u-channel inclusive scattering. 

In contrast, even were the on-shell approximation to be a good one, w-channel neutrino 
scattering does not allow one to extract the Gi at a single value of Q^. Since U"^ rather 
than is fixed, the value of varies as one integrates over the allowed region in {S,p) 
space. For example, for kinematics typical of the LSND experiment (see Fig. ^), varies 
over the range (see next Section) < \Q'^\ < 0.06 (GeV/c)^ as S and p vary over the 
allowed region V (allowing the neutrino scattering angle to vary over all possible values). 
Similarly, for kinematics typical of the BNL experiment p[, e = 1.3 GeV and, say, T/v = 500 
MeV, the corresponding range is < < 2.3 (GeV/c)^. In either case, the form factors 
associated with the struck nucleon contribute to the w-channel cross section over a range of 
Q^. In order to extract information about the form factors from this cross section, one is 
forced to adopt some parameterization for the Q^-dependence of the form factor and fit the 
parameters to the measured cross section. Furthermore, as alluded to above, in w-channel 
inclusive scattering the values of q, u and the neutrino scattering angle 6i, all vary when 
performing the integrations and thus no separation into isolated L, T, . . . hadronic responses 
is possible. 

All of these aspects of w-channel inclusive scattering stand in contrast with the t-channel 
situation discussed above in which a parameterization-independent form factor determina- 
tion is possible. Only in the case of elastic scattering from A = 1 targets are u- and t-channel 
processes equivalent. In the latter case, one has 

{K + PaY + + = 2m% + 2m^, , (36) 

where is the mass of the lepton and {K + Pa)^ = (e + mAr)^ — fe^^ m% — 2em7v in the 
lab frame. For scattering from a single-nucleon target, then, specifying f/^ is equivalent to 
specifying Q^. Consequently, ^-channel scattering can be used to perform a parametrization- 
independent form factor determination in this case. 

Finally, as discussed in Sec. and the Appendix, another issue to be confronted is the fact 
that in general the struck nucleon is initially off-shell and accordingly the relationship to on- 
shell single-nucleon form factors is not obvious. In the present work we use a generalization 



of the popular ccl prescription of de Forest |25| when modeling this type of behaviour. At 



least within the context of this approach we shall see that this model dependence is rather 
weak for the observables of interest (see Sec. |VlD . 
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IV. RFG PREDICTIONS FOR it-INCLUSIVE SCATTERING 



Although aheady considered in Ref. [Q, we wish to revisit in this Section the RFG 
predictions for the w-inclusive scattering partly to complement the findings of that work, 
which are rederived here via the alternative route of integrating in proper regions of the 
{S,p) plane, and partly to establish analytic expressions for some general features of the 
RFG M-inclusive cross section not presently available in the literature to our knowledge. 

As discussed briefiy in the previous Section, while in the t-channel it turns out to be 
possible to obtain a compact, simple expression for the RFG inclusive cross sections or, 
equivalently, for the longitudinal Rl and the transverse Rt response functions, this is not 
so in the w-channel. The reason is precisely the one mentioned at the end of the previous 
Section, namely, the Q^-dependence arising from the presence of the single-nucleon form 
factors. For the t-inclusive case, is fixed so that the nucleon form factors may be factored 
out of the integrals over S and p. By contrast in the w-channel varies over the integration 
region P; hence the nucleon's form factors cannot be brought out of the integrals except 
(and then only approximately) in special kinematic situations where the G'j(Q^) vary gently 
over V. 

To provide an appreciation for the behavior of in the {S,p) plane we display its 
variation in Fig. ^ for the same kinematical situation of Fig. ^, considering, as an example, 
the instance of Pn and k being parallel. The explicit dependence upon S and p of turns 
out to be 



Q' = {TN + Es + Sf (37) 

+ 



\kTPN 

where the upper (lower) sign corresponds to k and being parallel (anti-parallel). From 
Fig. ^, it clearly appears that varies rapidly with p but only mildly with £. 

Notwithstanding this feature we try to express analytically a few general properties of 
the RFG M-inclusive cross section. To this end, we recall |2y] that the RFG spectral function 
reads 

SRFG{p,S)=2e{kF-p) (38) 

X 5{8 - ^JPp + m^ + ^Jp'^ + m% 

where the factor 2 accounts for the spin degeneracy. Note that the above cannot be directly 
derived from the expression in Eq. (p9D, where the states are normalized according to (p4|). 
Instead, one requires the normalization 

{p\q) = 2Ep Q 6k,p , (39) 



n being the (large) volume enclosing the Fermi gas. Clearly, the support in Eq. (P5| ) in the 
{S,p) plane is nonzero only along the curve 



S^^^ = Jkl + m%- Jp^ + m% . (40) 
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It is then natural to expect the maximum of the w-inclusive cross section to correspond to 
the situation where T> encompasses the whole of the RFG spectral function. In the strict 
RFG model this requires 



k\ + — = ep — rriN = Tp > , (41) 

which insures that the piece of the boundaries of T> stemming from does not cut off a 
section of the curve defined by Eq. (0). 

For kinematic conditions under which the G't(Q^) vary mildly over the integration region 
T> and over variations in {uo,u), it follows from Eq. (|4l|) that the maximum of the cross 
section will occur for the following kinetic energy of the outgoing nucleon: 

^ l + m7v/(2A;) + A;/(2m^)sin2^jv ■ 

For purposes of comparison, we note that the same principle of maximum overlap between 
the support of the spectral function and T> leads to the following result for the maximum of 
the cross section: 



)2 I /r,™ — r^2 , ,2 



UJmax = \Qmax\/'^^N = [q -U;„J/2m^, (43) 



in the t-channel (see Ref. pO[). Actually Eq. (^) holds exactly for properly reduced 
RFG response functions where the single-nucleon form factor dependences are removed (see 
Ref. [^). If not so reduced, then at high momentum transfer the prediction in Eq. (^31) 



is appreciably altered by the single-nucleon physics. For the w-channel case, we find that 
the expression in Eq. (|4^), which provides the maximum for the cross section, is valid only 
for large k and small 9^, where it turns out that the single-nucleon physics has less of an 
impact. 

A novel feature of the RFG u-inclusive cross section, with respect to the t-inclusive one, 
is its unexpected vanishing in some range of T/y. To illustrate this feature, we recall that the 
M-inclusive cross section is fixed by three parameters: k, and 9^. Now for p~ = u — 6 > kp 
no overlap exists between T> and the curve of Eq. (|40|) : Hence the vanishing of the cross 
section. This situation clearly corresponds to a negative value of the intercept /„. However, 
when the intercept lu is positive, the cross section might also vanish, provided that 6 — u is 
larger than kp. One is thus led to consider the equation 



which admits two real, positive roots, namely 

ink'^ cos'^ 9f^ — a{2k — a\ 

2{k + rriN — a + kcos9M){k + ttin — a — kcos9N 



(44) 



(i_2) 2m„/c^ cos^ ^jv — Q;(2fc — a)(/c -|- — a) ± A;| COS 6'Ar|v^ . 



with 



A = 4m%k^ cos^ 9^ 

- a{2k-a){2mN-a){2k + 2mN-a) (46) 



14 



if, and only if, 



a(2m]\f — a)(2k — a)(2k — a + 2miq) 

^ w ■ 

where a = Es + kp and 

Es = mN- Jkl + ml ^ --^ (48) 

is the negative separation energy of the RFG (see Ref. [^]). 

It thus appears that a critical angle, 6*0, exists such that for 6^ < Oq the w-inclusive cross 
section vanishes in the range T^y"* < < Tj^^ . The asymptotic value of Oq for large lepton 
momenta k is 52° at kp = 225 MeV/c and exists only for k > kp. Note that when k = kp 
then On = 0^ and the two roots in Eq. (^) coincide, namely 

= = -Es . (49) 



For k > kp the maximum distance between the roots in Eq. (the maximum range over 
which the cross section vanishes) occurs for 9n = 0. Then it gradually decreases as 9n is 
increased from 0° to Oq. 

A final consideration relates to the traditional handling of Pauli blocking in the RFG. 
In the t-channel, this blocking gives rise to the experimentally unsupported linear energy 
behaviour of the cross section at low u when q < 2k p. Consequently, one does not expect 
the RFG cross section to be credible until q > 2kp. In the u-channel, Pauli correlations 
may also be incorporated, although one may debate the appropriate method for treating 
them. At the crudest level of approximation, for example, one may model these correlations 
simply by requiring that the cross section vanish for < kp (T/v < Tp). Importantly for 
comparisons with the LSND measurements, an experimental cut is made at T/v = 60 MeV 
(Pat ~ 340 MeV), lessening the impact of such modeling. 

Let us now conclude this Section by presenting some typical w-inclusive neutrino cross 
sections using the RFG. For this purpose we first recall that the RFG is a symmetric system 
with equal number of protons and neutrons (Z=N=A/2), all of them on the mass-shell. 
Then we insert the spectral function of Eq. (^) into Eq. (|35|) , yielding 

(50) 



dVLNdEN p(27r)n6 u KMy 

.Srpg{p,S) 



where p = 2fc|./(37r^) is the Fermi gas density and N" indicates either the proton (Z) or 
the neutron (N) number. Furthermore, when applied to neutrino scattering, the coupling 



to mass ratio which enters is g/{4:Cos6wMw) = \/Gf/{2V2), Gp being the Fermi constant 
(we have assumed the p parameter to be equal to unity for simplicity), and |Q 

L^.w^'ip, pn) = 2[Vvv + VvA + Vaa] , (51) 
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where explicit expressions both for the on- and off-shell quantities in Eq. (|5T|) are given in 
terms of the weak neutral current form factors in the Appendix. 

Since our focus in this work is on ^^C, we take kp = 225 MeV/c. The results of our 
calculations are displayed in Fig. where we again consider kinematics typical of the 
LSND experiment (e = 200 MeV) and, for purposes of illustration, we explore two different 
orientations between the incoming neutrino and the outgoing nucleon (either proton or 
neutron), namely 6n = 20° and 60°. We observe the following: 

i) the cross section decreases with the angle 9n for low values of k, at least for not too 
large values of 6' at; 

ii) the neutron cross section is larger than the proton cross section. 

In Fig. 1^ the neutrino cross sections for the RFG (but with only outgoing protons) are 
calculated using kinematics corresponding to the average energy of the Brookhaven neutrinos 
(e = 1.3 GeV). In this case, the experimental cut is made at T/v = 200 MeV and so Pauli 
blocking plays a minor role. In contrast to the situation for low-energy neutrinos, the cross 
sections increase with 9n- In addition, the previously noted possibility of a vanishing u- 
inclusive cross section for certain values of 6n appears in Fig. |^: A sizable gap in the cross 
section occurs for 9n = 20°, but not for 9^ = 60°, since the value of the limiting angle in 
this case turns out to be 9o = 25.7°. 



V. HM PREDICTIONS FOR u-INCLUSIVE SCATTERING 



We next turn to the hybrid harmonic oscillator shell model 0] and its extensions to 
inclusion of a spreading width. While these choices clearly do not exhaust the list of possi- 
bilities — in particular in this work we do not deal with the final-state dynamics issues that 
were treated in Ref. |T^ — we are able to set a scale for the theoretical uncertainty in the 
extracted value of GyT as discussed in the next Section. 



The HM developed in Ref. |^0[ differs from the RFG in that it has the states below 
the Fermi surface confined, therefore requiring the struck nucleon to be off-shell (see the 
Appendix). In the PWIA the states above the Fermi surface are taken to be plane waves. 
Accordingly in this model one accounts for the confinement of the nucleons in the initial 
state, which is of course important for a realistic description of low-energy neutrino-nucleus 
inelastic scattering, however at the expense of losing the Lorentz covariance characteristic 
of the RFG. In other words the hybrid model is only semi-relativistic. 

In principle, this confinement is generated by a Hartree-Fock well; in practice, as a simple, 
tractable model we use a harmonic oscillator as the confining potential, namely 

Vir) = ^rriMOoy - V , (52) 

V being a positive quantity fixed in such a way to reproduce the experimental separation 
energies. Eg = 15.96 (18.72) MeV for protons (neutrons). For the oscillator frequency, we 
use ujq = 41/y4^/^ MeV. The HM spectral function is then easily found to be 
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Nmax 

S''^'{p,S)= J2 S[S-{NMAx-N)uo + AE]nM 

N=0 

Nmax o 
Af=0 ^ 



(53) 



where the shift in the excitation energy £ stemming from the binding of the nucleons is 
denoted by AE. In the above, the momentum distribution for a given shell identified by the 
quantum number N = 2n + C. given by 

9/ -I- 1 

nM=2Y^^\vM\\ (54) 

where <fne{p) are the harmonic oscillator radial wave functions in momentum space. For 
simplicity (and in good accord, for example, with mean-field descriptions of light N = Z 
nuclei) we assume the momentum distributions of protons and neutrons in a given shell to 
be the same. In contrast with the RFG spectral function discussed above one finds that the 
HM spectral function in Eq. (|53|) is nonzero on a set of parallel lines in the {S,p) plane, one 
for each shell (e. g., two lines for ^^C), at variance with the case of the RFG, whose spectral 
function is nonzero on the single curve given by Eq. (|3). As mentioned in the Introduction, 
in the present work when considering the HM we employ only a harmonic oscillator basis 
with its overly strong confinement to provide the strongest contrast with the RFG model — 
a sort of "worst-case" scenario. 

As one straightforward extension of the above formalism where specific classes of corre- 
lations effects may be modeled we also allow the single-particle states to acquire a spreading 
width. This extension is implemented by adding a complex self-energy to the single-particle 
energies, thus defining the fermion propagator as 

GS'>(p, E) = T , , ""ff , , , (55) 



N=0 



E -[eN + S(eAr) - ep] 



where = rn^ + (A^ + ^)ujo — V represents the energy of the nucleon in the A^th shell. 
The self-energy S, of course, arises from second- (or higher-) order insertions on the fermion 
propagation lines: Here, we adopt the point of view of parametrizing it in terms of the 
following function 

S(e) = A(e) - , (56) 

with 

7W = 2«;^;^^(-^) , (57) 
e -|- to t -l- ti 

whereas A(e) is obtained from 7(e) via a subtracted dispersion relation. We follow 



Ref. |32,3l in setting a = 10.75 MeV, eo = 18 MeV and ei = 110 MeV. The spectral 
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function is proportional to the imaginary part of the hole propagator, S\j^.i = ImGjjj^ /ti: 
Hence, in principle, the spreading width extends the support of the spectral function to the 
whole first quadrant of the {S,p) plane; in practice, S^^^ remains concentrated around the 
lines where the spectral function for a pure harmonic oscillator is non- vanishing. 

We turn now to an examination of some of the model dependences that we obtain. In 
Figs. 1^ and ^, we display the neutrino cross sections for the HM and its extended version with 
a spreading width for the same kinematical conditions explored above using the RFC One 
sees from the figures that, although the basic features of the RFC cross sections studied 
in Section |V] do not appear to be much altered by the nucleons' confinement, yet some 
differences do show up. In particular, the peaks of the cross sections in the HM appear 
shifted and the phenomenon of the vanishing of the cross section, which is characteristic of 
the RFG, is still present in the hybrid model, although in a narrower range of energies. Note 
also in comparing panels b) and c) that the effect of the spreading width, and therefore this 
particular class of correlation effects, is negligible and is henceforth ignored. 

While the nuclear modeling used in these comparisons of RFG and HM results can ob- 
viously be extended in many directions, including the use of alternative mean-field bound 
wave functions and incorporation of final-state interaction effects, the theoretical uncertain- 
ties are likely reasonably characterized by the present "extreme" situations. Each model has 
its own merits. In particular, the RFG maintains covariance, which becomes more relevant 
for the BNL kinematics, whereas the HM brings in the confinement of the (off-shell) struck 
nucleon, albeit at the expense of covariance. Since the former involves only plane- wave nucle- 
ons while the latter has rather strongly confined HO wave functions for the struck nucleons, 
most other choices might be expected to fall somewhere between the results presented here. 
If the span of predictions given in the figures for the double-differential cross section is taken 
at face value, then substantial nuclear model dependences are anticipated and any hope of 
using absolute cross sections for the purpose of learning about axial-vector strangeness in 
the nucleon is correspondingly low. As we discuss in the next section, however, there does 
appear to be a strategy whereby one can mitigate this model-dependence in the extraction 
of As from neutrino observables. 



We now address the use of neutrino scattering to probing nucleon strangeness. To this 



VI. STRANGE AXIAL- VECTOR FORM FACTOR 




where 



Gf = {l/2){D + F)GUr) 
Gf = {l/2V3){3F-D)GUr) 



(59a) 
(59b) 
(59c) 
(59d) 



D + F = 1.257 
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are SU(3) octet axial-vector form factors, D and F are the associated SU(3) reduced 
matrix elements, g'^^ = As is the strange quark contribution to the nucleon's axial charge, 
and r = |Qp/4m^. The coefficients (^^^ are determined by the axial- vector coupling of 
the to the quarks, and take on, at tree level in the Standard Model, the values —2, 0, 
and 1 for a being T = 1, T = 0, and 0, respectively. For neutrino scattering, electroweak 
radiative corrections are of (9(a/47r) [^,0, so for purposes of this analysis, we may employ 
the tree-level values. Under this approximation, any isoscalar component of the nucleon's 
axial- vector form factor arises solely from the strange quark term. 

The quantity G^{t) is a dipole parameterization for the Q^-dependence of the axial- 
vector form factors. In principle, one has no rigorous justification for choosing a dipole form 
and for assuming that the different form factors in Eq. ( p^D display the same Q^-behavior. 
In the case of high-energy neutrino scattering, where values for \Q'^\ on the order of one 
(GeV/c)^ are achieved, these assumptions can have rather drastic consequences for one's 
determination of As. As indicated in Refs. 0J^, the value of As is strongly correlated with 
the dipole mass parameter. Ma- As far as a high-energy neutrino scattering determination 
of As is concerned, the impact of choosing a different parameterization or of allowing the 
octet and singlet form factors to display different Q^-behavior is unknown. At kinematics 
relevant to the LSND experiment, the issue of non-leading Q^-dependence is much less 
serious. In this case, one may view the dipole form factors as a way to include the first 
derivative-term in a Taylor expansion of the form factor (with respect to Q^). Moreover, as 
shown in Refs. 0J^, the correlation between As and Ma is negligible. Thus, we may safely 
employ the dipole parameterization for purposes of analyzing low-energy scattering without 
introducing significant uncertainty. 

To illustrate the relationship between As and the choice of nuclear model, we first plot 
in Fig. ^ the differential cross section as a function of outgoing nucleon kinetic energy for 
two different values of As: and —0.2. The latter value corresponds roughly to the largest 
value (in magnitude) of As derived from the deep inelastic measurements. The results for 
kinematics typical of the LSND experiment are shown in Figs. Pa,b while those typical 
of the BNL measurements are given in Fig. |^c. We note that in the case of the LSND 
experiment, for which the target is a CH2 based scintillator and for which the incident 
neutrinos are produced by the decay of pions in fiight, a cut on T/v of > 60 MeV guarantees 
that only nucleons knocked out of carbon nuclei are detected ||12|. For these energies, one 
sees from Figs. Ha,b that the nuclear model-dependence is sufficiently large to prevent an 
unambiguous extraction of As. For example, in the case of the proton knockout cross section, 
the difference between the RFG prediction for As = and the HM prediction for As = —0.2 
is comparable to or smaller than the difference between the RFG predictions with As = 
and As = —0.2. This feature does not appear to persist at higher-energies (Fig. ^) where 
the model-dependence for a given value of As is significantly smaller than the dependence 
on the value of As. 

The authors of Ref. |]12| proposed that instead of analyzing the proton and neutron cross 
sections separately, one ought to consider the ratio of total proton to neutron yields. In 
contrast to the situation with the individual cross sections, the yield ratio is insensitive to 
one of the primary experimental uncertainties - the normalization of the incident neutrino 



fiux. Moreover, it was shown in Ref. that the yield ratio is less sensitive to final- 
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state interactions than are the individual cross sections. Thus, the use of the p/n ratio 
would appear to minimize the impact of experimental and theoretical uncertainties on the 
extraction of As from QE neutrino scattering data. In what follows, we illustrate the degree 
to which the use of this ratio can also reduce one's sensitivity to the choice of nuclear model. 



We work with the ratio of cross sections rather than yields, assuming as in Ref. |T2[ that 
the incident fluxes are peaked at some energy e and that working with the cross sections 
at e = e is sufficient. To that end, we first show in Fig. ^ the ratio of differential proton 
knockout and neutron knockout cross sections as a function of T^r. From these curves, we 
deduce that the spread due to the choice of nuclear model for a given value of As and Tjv is 
roughly 20% or less than the change in the ratio obtained by varying As from zero to —0.2. 

In Fig. |lO| we give the total cross sections, integrated over outgoing nucleon energy (with 
Tjv > 60 MeV in the case of low-energy neutrinos and T/v > 200 MeV for high-energy 
neutrinos) as functions of As. In the case of the LSND experiment, it is the total proton 
and neutron yields that will be obtained, making these curves more directly relevant to the 
interpretation of the experiment. From Fig. p!0|a,b, one sees that the model-sensitivity of 
one's extracted value of As is non-trivial in the case of low-energy neutrino scattering. For a 
given value of the total cross section, the two would yield values of As differing by about 0.5. 
Thus, we would take a reasonable nuclear theory error bar on the value of As extracted from 
the single cross section to be 5nuc(^'5) = ±0.25. This uncertainty is as large in magnitude 
as the value of As determined from deep inelastic scattering. For higher energy scattering, 
however, the corresponding nuclear model uncertainty is much smaller: (5nuc(As) = ±0.015 
(see Fig. |l( 



Finally, in Fig. [11] we give the ratio Ry of proton to neutron total cross sections for low- 
energy scattering. In this case, one finds that the ratio is significantly less sensitive to the 
choice of nuclear model than in the case of the actual cross sections themselves. Indeed, we 



would deduce a nuclear theory uncertainty in As of ~ ±0.015 from the curves in Fig. |r^. This 
uncertainty represents an order of magnitude improvement from the uncertainty associated 
with the use of the individual cross sections. 

One other issue involving model dependence was also briefiy examined, namely, that of 
the on- versus off-shell current descriptions employed. We evaluated the HM results using 
the on-shell single-nucleon current and found that Rj^ in that case is lower than the RFG 



curve by about the amount it is higher for the off-shell results shown in Fig. ITT|. Of course. 



such a calculation is not strictly correct and should be used only as a rough indication of the 
sensitivity to the prescription for the current, since some off-shellness must be used for the 
HM as demanded by the kinematics of the reaction. In fact, the rather weak dependence 
found is reassuring and suggests that, at least for the specific kinematical choices that are 
relevant here, the issue of what particular single-nucleon current to employ is likely even less 
important than the nuclear model dependence. To be as conservative as possible, we will 
multiply by two the uncertainty deduced from Fig. [Tl| and take 5nuc(^'S) = ±0.03. 

The reason for this reduction in nuclear model sensitivity when using R^ is straightfor- 
ward to understand. At the kinematics typical of the LSND experiment, the range of 
accessed as one integrates over the allowed region V^On) is small. Consequently, the nucleon 
form factors and kinematic factors vary gently and, to a first approximation, may be fac- 
tored out of the integral in the expression of Eq. (|35D. One is then left with a product of the 
squared single-nucleon invariant amplitude and an integral of the one-body spectral func- 
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tion. Differences between model predictions would then arise from different distributions of 
strength in the spectral function over the allowed domain in the {S,p) plane. In the ratio 
of proton and neutron cross sections, the integrals of the spectral functions should cancel 
assuming the proton and neutron S{p, £) are identical. In this case, the ratio is independent 
of nuclear model and is determined only by the single-nucleon form factors. Any model de- 
pendence which appears in the ratio would then arise from one of the following sources: (i) 
the variation in form factors and kinematic factors over the integration region V{6n), so that 
the factorization into a product of single-nucleon and many-body physics is not exact; (ii) 
differences between proton and neutron spectral functions; (iii) contributions going beyond 
the PWIA. In the case of low-energy neutrino scattering from carbon for T/v > 60 MeV, one 
expects the effects (i) and (ii) to be small. Indeed, our results indicate that effects of type 
(i) generate roughly a 10% variation in the ratio for a given value of As, which translates 
into the nuclear theory uncertainty quoted above. 

In arriving at our estimate of the nuclear theory uncertainty, we make no pretense of 
having performed definitive, state-of-the-art nuclear model calculations of the kind reported 
in Refs. |T2| , p!3[ . Rather, our goal was to give a reasonable upper bound on the expected 
spread in model calculations. We anticipate that any series of more sophisticated model 
calculations would yield extractions of As from differing by no more than our value for 
^nuc(^'S). From this standpoint, it is instructive to compare our RFG and HM results for 



this ratio with the mean field calculation reported in Ref. |T^. For a given value of As, the 
mean field results and RFG results for are almost identical, while the HM results differ by 
less then 10%. This agreement holds in spite of much larger differences which appear when 
one compares the mean field, RFG, and HM predictions for the individual differential cross 
sections and for the ratio of differential cross sections as a function of T/y. We further believe 
that this robust nature of R^, will persist when one compares different model calculations 
which not only reproduce experimental results for the inclusive t-channel QE responses but 
also include FSI's in the w-channel case. It was demonstrated in Ref. that the inclusion 
of FSI's in the mean field approach can increase the predicted value of R^ by more than 
10% over the mean field and RFG predictions. This increase would imply a shift in the 
extracted value of As by more than 0.03. Thus, it appears that any realistic model which 
is used in this extraction must include FSI's. Nevertheless, we would argue that the spread 
in extracted values of As corresponding to the use of different nuclear models which include 
FSI's would remain smaller than 5nuc(As). Indeed, an important check on our first estimate 
of ^nuc(^'S) would be to compare model calculations which incorporate FSI's as well as to 
explore the impacts of correlations and charge symmetry breaking. 

Finally, we compare the nuclear theory uncertainty in the neutrino scattering deter- 
mination of As with the theory uncertainties associated with the extraction of As from 
measurements of the gi{x,Q'^) deep inelastic structure function. We first remind the reader 
that one does not necessarily expect the values of As obtained from these two different pro- 
cesses to agree. It has long been known, but perhaps not widely appreciated in the nuclear 
physics community, that As is a renormalization-scale dependent quantity. In deep inelastic 



scattering, this scale is typically taken to be the same as the mean y\Q'^\ of the reaction 

(between ^/2 and v^TT GeV/c), whereas the scale appropriate to the quasielastic neutrino 
value is somewhere below the charm quark mass. The latter corresponds to the scale below 
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which one includes only the three lightest quarks in an effective axial- vector NC [|3J]. It is 
often assumed that the QCD evolution of As between the two scales is rather gentle, based 
on leading-order perturbative calculations, yet it is conceivable that non-perturbative effects 



could invalidate this assumption |3^. At present, one can make no definitive statements 
regarding the evolution of As between these two scales. 

A separate issue pertaining to the use of polarized structure functions to determine As 
is the use of SU(3) symmetry. In the standard operator product analysis of the (71-sum, one 
employs an SU(3) parameterization of the octet axial- vector matrix elements {N\A^^\0)\N) 
(a = 3, 8 refers to the respective components of the octet) in order to derive a value for 
As or AS. The latter denotes the nucleon's singlet axial charge or, in the quark-parton 
framework, the total light quark {u, d, s) contribution to the nucleon spin. The value of As 
determined in this fashion is quite sensitive to the quantity 

-|(iv|^(8)(o)|iV) = l(3F-D) . (60) 

In an SU(3)-symmetric fit to hyperon semileptonic decays, the combination of reduced ma- 
trix elements 3F — D takes on a value of about 0.6 [|23[. Because it involves a cancellation 
between two quantities, this number is rather sensitive to uncertainties in the fit as well as 
to corrections arising from SU(3)-breaking in the octet matrix elements. Several analyses 



have been carried out recently ||22|j2^,|36|j33] in which the hyperon decays were re-fit with al- 
lowance for SU(3)-breaking. While different schemes for the incorporation of SU(3)-breaking 
were employed in each case, a general trend does emerge: SU(3)-breaking may reduce the 
matrix element in Eq. (pOD by 50% or more compared to its SU(3)-symmetric value, with 
an uncertainty of comparable magnitude. Such a reduction would imply a value of As ~ 
from deep inelastic data, in contrast to the current average of the measurements As ~ —0.1. 
The authors of these studies caution that precise numerical value for the SU(3)-breaking 
correction to the quantity in Eq. ( pDj ) is not highly reliable, due to the nature of the afore- 
mentioned cancellation. To be conservative, then, we take the SU(3)-breaking uncertainty in 
this matrix element to be 50% of its SU(3)-symmetric value. The corresponding uncertainty 
in the strange-quark axial charge is 5]2)jg(As) ~ ±0.1, a value having the same magnitude 
as the present average for As under the assumption of good SU(3) symmetry in the hj^eron 
semi-leptonic decays. The scale of this uncertainty is consistent with the scale of the error 



obtained with the heavy-baryon chiral perturbation theory analysis of Ref. p2| and with 
the range of estimates for SU(3)-breaking corrections given in Refs. [^,^,^. The presence 
of this SU(3)-uncertainty weakens the standard conclusion drawn from inclusive, polarized 
deep inelastic measurements that the strange quarks are polarized oppositely to the direc- 
tion of the nucleon's spin and that the magnitude of the ss contribution is about one-third 
the total quark contribution, AS. 

By way of comparison, we note that ^nucl^-s) is about one-third as large as 5£)jg(As). 
Moreover, the analysis of the QE neutrino cross sections does not suffer from the kind of 
SU(3)-breaking corrections and uncertainties which enter the interpretation of the gi-sum 
results. Although the same problematic combination of SU(3) reduced matrix elements 
3F — D enters the decomposition of the nucleon's axial- vector NC form factor (see Eqs. 
59]) , its contribution is suppressed by the coupling The latter is identically zero at 



tree level in the Standard Model and is on the order of 0.01 when one-loop electroweak 
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corrections are included. From this standpoint, then, it appears that the ratio Ry provides 
a theoretically cleaner window on the strange-quark contribution to the nucleon's spin than 
does inclusive, polarized deep inelastic scattering. The primary limitation in the former case 
appears to be experimental error. 

VII. CONCLUSIONS 

The study of low- and intermediate-energy semi-leptonic processes has played an impor- 
tant role in helping uncover the nature of fundamental interactions. In the case of QCD, 
semi-leptonic scattering is poised to illuminate the way the strong interaction is realized in 
the structure of the nucleon. As with the use of semi-leptonic scattering to test the Stan- 
dard Model and its possible extensions, the efficacy of this probe of the nucleon's structure 
is limited by the reliability with which one can compute theoretically, or determine experi- 
mentally, the other hadron and nuclear structure contributions to the relevant observablcs. 
It is important, then, that where hadron structure or many-body theory is brought to bear 
on the interpretation of such observables, an attempt be made to quantify the theoretical 
uncertainty. 

In the present work, we have attempted to provide a theoretical error bar for the deter- 
mination of As from quasielastic neutrino NC scattering. The nucleon's strange quark axial 
charge is of interest since it has the quark-parton interpretation as giving the strange-quark 
contribution to the spin of a polarized nucleon. A large value for this quantity would imply 
that the non- valence quarks play a more central role in the low-energy characteristics of the 
nucleon than implied by the highly-successful and intuitively satisfying quark model. To 
the extent that one may extract the nucleon's NC axial- vector form factor from quasielastic 
neutrino scattering at low momentum transfer, one has a direct probe of As. In the ideal 
situation, one would have in hand sufficient information from quasielastic (e, e'N) measure- 
ments to perform this extraction without heavy reliance on theoretical calculations of the 
QE response. At present, sufficient (e, e'N) data are lacking, rendering the use of nuclear 
models unavoidable. 

We have tried to argue that were one to rely on the individual {v^N)v' cross sections 
measured by LSND, the nuclear theory error bar on As would be sufficiently large to render 
the results inconclusive. Fortunately, the ratio of total proton to neutron yields {R^,) appears 
to be much less sensitive to the choice of nuclear model, reducing by nearly one order of 
magnitude the nuclear theory uncertainty in As. Based on the study of some simple nuclear 
models with quite distinct assumptions regarding the many-body dynamics, we estimate 
the nuclear theory error to be 5nuc(As) = ±0.03. We anticipate that the spread in values 
for As obtained from R^j and more sophisticated nuclear models will be smaller than this 
value for the uncertainty. We further note that 5nuc(As) is about a factor of three smaller 
in magnitude than 5£)jg(As), where the latter is derived from recent analyses of SU(3)- 
breaking in the octet of baryon axial- vector matrix elements. In contrast, the interpretation 
of QE neutrino cross sections is relatively free from such SU(3)-uncertainties. 
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APPENDIX: 



The quantities appearing in Eq. ( |^ ) are given by the following expressions: 

Vvv = '2[Vn + Vi2 + V22] 
Vva = 2[Vai + Va2] 
Vaa = 2A 

using the nomenclature of Ref. [Q. On-shell one has from that reference 
Vii = AF^{P ■ K Pn ■ K' + Pn ■ K P ■ K' - m^K ■ K') 



(Al) 
(A2) 
(A3) 



(A4) 



12 



-AF1F2K ■ K' {Pn ~P)-{K - K') 



(A5) 
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2F| 



m 



K ■ K' {P ■ K Pn ■ K + P ■ K' Pn ■ K' 



N 



+ m%K ■ K') 



A = AG\{P ■ K Pn ■ K' + Pn ■ K P ■ K' + m%K ■ K') 



(A6) 
(A7) 



Al 



^GaFi{P ■ K Pn ■ K' - Pn ■ K P ■ K') 



(A8) 



Va2 = AGa^K -K' {P + Pn) ■ (K + K') 



(A9) 



where P and Pn are the initial and final four-momenta of the nucleon, K and K' the initial 
and final four-momenta of the neutrino, Ga is the axial-vector form factor of the nucleon 
whereas Fi and F2 are the Dirac and Pauli weak NC form factors of the nucleon, related to 
the Sachs NC form factors through the expressions Ge = Fi— TF2 and Gm = A + -^2, with 
= IQ^I/^'^AT- For the nucleon's form factors we have used the Galster parameterization 
(see Refs. [p8|pT| ). These can be combined as in Eqs. ( [A1| -[A^) to obtain the following for 
the on-shell versions of the quantities in Eq. (pl|): 



TAon 

yvv 



^VA 



^AA 




16r(G| 



tGmGa 



16r(l + r) 



g\ 



(AlO) 
(All) 
(A12) 
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P,Y = K + K' and 
= 0. Furthermore, 



we 



Here we have used the following four-vectors: Q = K — K' = P/v - 
Z = P/v + P, where then Q ■ Y = and for on- shell nucleons Q ■ Z 
define as usual W2 = [G| + tGIj]/ {1 + t). 

Off-shell we use a generalization of the ccl prescription introduced by de Forest ||25 
for electron scattering; in our case we now have analogous quantities involving both vector 
and axial- vector neutral currents. Defining, as usual, = Pj^ — = {E,p), letting 
pi' = P^ - Qi' = (E,p), with E = + as in Ref. ^ 

Z'^ = Pj^ + P^, so that Q ■ Y = Q ■ Z = 0, we obtain 



25 1 , and in addition defining 



yoff 
VVV 



+ 



yoff 



yoff 
^AA 







Gl, + AmlQ^G'l 


(F-z)'-(g 






1 [{q ■z){q- y) 




(q-q) {y-z) 


GmGa 


(Q.Qf+(Y 


-\2 
^) 




(g-z)'-(g- 


\2 
^) - 


- Am\Q^\ G\. 



(A13) 

(AM) 

(A15) 



In Eq. ( [A13| ) we have introduced the following additional single-nucleon form factors: 



G' = Pi - rP2 



1 



G'J + tG\,i 



(A16) 
(A17) 



in close analogy with the above on-shell quantities (see also Ref. ^8[). The results given 
here for the half-off-shell currents revert to the on-shell results in Eqs. ( |A1(]| - |AT2D when P 
and P become the same, i.e., when the kinematics force the nucleon on-shell. 
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FIGURES 



FIG. 1. Feynman diagram for exclusive lepton scattering in the one boson exchange approxi- 
mation. 

FIG. 2. The domain T> over which one integrates for the f-inclusive cross section at typical 
values of momentum and energy transfer. Note the more complex structure of the boundary when 
the intercept It given by Eq. || is positive. 

FIG. 3. The same as Fig. 2 for the li-inclusive cross section: a) LSND kinematics (e = 200 
MeV, Tn = 60 MeV); b) BNL kinematics (e = 1.3 GeV, Tn = 60 MeV). The boundaries shown 
all involve £~ except for line BF which involves £^ . Note that when the neutrino and outgoing 
nucleon momenta are antiparallel quite remote regions of the (<f,p) plane are explored. 

FIG. 4. Plane- wave impulse approximation version of the diagram in Fig. 1. 

FIG. 5. Variation of the squared four-momentum transfer in the {£,p) plane for the kinemat- 
ical conditions of the LSND experiment. Note the rather mild dependence on £ and the strong 
dependence on p. 

FIG. 6. The n-inclusive double differential cross section for protons (solid) and neutrons 
(dashed) for LSND kinematics and two different orientations of k and pN- (a) RFG model, (b) 
HM with harmonic oscillator wave functions and (c) as for (b), but including a spreading width, or 
equivalently, a complex nucleon self-energy. Note the negligible effect of the last. Note also that, 
although we display the results over the whole allowed range of T/v, in the LSND experiment only 
nucleons with T/v > 60 MeV are actually detected. 

FIG. 7. The same curves as in Fig. 6, but now for BNL kinematics. Note the double peak 
behaviour of the cross section at ^at = 20^ (see text for discussion). Note also that, although we 
display the results over the whole allowed range of T/v, in the BNL experiment only nucleons with 
T/v > 200 MeV are actually detected. 

FIG. 8. (a) The effect of axial-vector strangeness in the RFG and HM proton ejection cross 
sections for LSND kinematics integrated over angle for two values of g\ . Note the significant 
increase of the cross section induced by strangeness, (b) The same as (a), but now for neutrons. 
Note the significant reduction of the cross section induced by strangeness, (c) The same as (a), 
but now for BNL kinematics. 

FIG. 9. Ratio between the angle-integrated cross sections for proton and neutron ejection at 
LSND kinematics. Both the RFG (dashed) and HM (solid) results are displayed. 
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FIG. 10. Total cross section (integrated over angle and energy) versus the strangeness content 
of the nucleon for outgoing protons (a) and neutrons (b) at LSND kinematics and (c) for protons 
at BNL kinematics. Dashed lines: RFG; solid: HM. 

FIG. 11. The ratio of the proton and neutron total cross sections for LSND kinematics: 
Dashed line: RFG; solid: HM. 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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Fig. 5 
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Fig. 9 
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